Primordial Gravity Waves and Weak Lensing 
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Inflation produces a primordial spectrum of gravity waves in addition to the density perturbations 
which seed structure formation. We compute the signature of these gravity waves in the large scale 
shear field. In particular, the shear can be divided into a gradient mode (G or E) and a curl mode 
(C or B). The former is produced by both density perturbations and gravity waves, while the latter 
is produced only by gravity waves, so the observations of a non-zero curl mode could be seen as 
evidence for inflation. We find that the expected signal from inflation is small, peaking on the 
largest scales at 1(1 + l)Ci/2ir < 10~ n at I — 2 and falling rapidly there after. Even for an all-sky 
deep survey, this signal would be below noise at all multipoles. Part of the reason for the smallness 
of the signal is a cancellation on large scales of the standard line-of-sight effect and the effect of 
"metric shear." 



The theory of inflation was proposed over twenty years 
ago[l]. For the first years after its discovery/invention, 
cosmologists worked out some of its cosmological impli- 
cations and particle physics implementations. While this 
work is still going on, the most exciting development 
in the last several years has been the confirmation of 
some of inflation's basic predictions. It now appears that 
the universe is flat, the most robust prediction of infla- 
tion [2]. Further, observations of large scale structure and 
anisotropics in the cosmic microwave background (CMB) 
strongly support the idea that small, adiabatic, nearly 
scale-invariant fluctuations were present in the early uni- 
verse. Inflation predicts the existence of precisely this 
class of perturbations [3] . 

The case for inflation is not airtight though. Infla- 
tion explains why the universe appears flat today even 
if the curvature is not exactly equal to zero. But it 
is possible that the curvature really is exactly zero, so 
that inflation is not needed to account for the appar- 
ent flatness. Similarly, the adiabatic, scale-free pertur- 
bation spectrum predicted by inflation might have been 
laid down by some other mechanism. Indeed, it is some- 
what of an embarrassment for proponents of inflation 
that this spectrum is called the "Harrison-Zel'dovich- 
Peebles" spectrum, named after the three eminent cos- 
mologists who first proposed it, long before inflation had 
been suggested. 

Are there any signatures unique to inflation? Just as 
the primordial density perturbations were produced by 
quantum fluctuations during inflation, primordial grav- 
ity waves were produced by quantum fluctuations of 
the metric during the inflationary epoch [4] . If we were 
to detect these gravity waves, the case for inflation 
would be strengthened considerably. After the detection 
of anisotropics in the CMB by the COBE satellite [5] , 
there was hope that the gravity waves would eventually 
be extracted from an accurate measurement of the full 
anisotropy spectrum [6] . Recent work has shown though 
that other cosmological parameters, especially late reion- 
ization, can mimic the effects of gravity waves in the 



anisotropy spectrum [7], making detection more difficult. 

In 1996, two groups [8, 9] showed that gravity waves 
produce a polarization pattern in the CMB that cannot 
be caused by scalar (density) perturbations to first order 
in the amplitude. The polarization field can be decom- 
posed into two modes. One is technically a scalar pattern 
on the sky, just as a gradient is, and the names given to 
such modes have variously been "scalar" , "G" (for gra- 
dient) , or "E" (for electric since an electric field is a 
gradient of a scalar potential). Another orthogonal set 
of modes are pseudoscalar patterns on the sky, and can 
be obtained by rotating the polarization at each direc- 
tion in a scalar pattern by 45°. Such modes change sign 
under parity like a curl does, and names given to them 
have variously been "pseudoscalar", "C" (for curl), and 
"B" (for magnetic since a magnetic field is also odd un- 
der parity). The polarization pattern produced by a sin- 
gle sinusoidal 1 density perturbation must be rotationally 
symmetric about the wavenumber k and also symmetric 
under reflections in the plane perpendicular to k. Such 
a pattern cannot contain any C (=B) modes. In lin- 
ear theory the polarization pattern produced by density 
perturbations is just the sum of the pattern from individ- 
ual sine waves and hence we derive the general rule that 
scalar perturbations produce no C modes 2 . This sym- 
metry argument does not apply to gravity waves which 
do not have the same reflection symmetry because each 
gravity wave is polarized in a particular direction. Grav- 
ity waves do produce C mode polarization patterns in 
the CMB. Thus detection of C modes can provide fairly 
unambiguous evidence for gravity waves, thereby verify- 
ing a unique prediction of inflation. It will take quite a 
while to reach the sensitivity needed to make this test a 



1 or pseudo-sinusoidal in the case of an FRW cosmology with non- 
zero curvature 

2 Nonlinearities can lead to C modes even from scalar inhomo- 
geneities but will be most important on small angular scales. 
The gravity wave signal discussed here is on large angular scales. 
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powerful probe of inflation [10]. Along the way, there are 
many systematic effects which might prove to be spoilers, 
prominent among them the possibility that foregrounds 
are polarized. 

The pattern of shearing of images due to gravitational, 
or any other kind of, lensing is described by a symmetric 
traceless tensor field on the sky, just as is the polarization 
pattern, and like polarization can be decomposed into G 
and C modes. The gravitational field causing the lensing 
can be scalar (caused by density inhomogeneities), vec- 
tor (caused by vorticity), or tensor (gravity waves). The 
symmetry argument stated above tells us that C mode 
shear pattern can be produced by gravity waves and not 
by density inhomogeneities [11]. As with CMB polariza- 
tion, we expect the dominant shear to be G modes pro- 
duced by density inhomogeneities, but we can hope to 
detect gravity waves through C modes in the shear pat- 
tern. This raises the question of whether the primordial 
gravity waves produced by inflation might be detectable 
by measuring the C mode of shear. Here, we explore this 
question. 

One measures the shear field by correlating the shapes 
of distant galaxies, in particular the ellipticity of galactic 
light distribution on the sky [12], and as with polarization 
we can expect various other effects to mimic C mode 
shear from gravity waves in such a measurement. We will 
not address the issue of "foregrounds" in this paper, but 
as we will show (see also Ref. [13]) the amplitude of shear 
produced by gravity waves is so small that the expected 
signal is below even the statistical noise expected in the 
largest experiment imaginable. 

Gravitational lensing is caused by the deflection of light 
trajectories due to metric perturbations. The metric 
of a flat Friedman-Robertson- Walker cosmology in syn- 
chronous coordinates 3 is given by 
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where a (r/) is the scale factor of the universe; rj is the 
conformal time variable; and the 3x3 matrices I and H 
are respectively the identity matrix and the metric per- 
turbation which is symmetric H = H T . The deflection 
of light is described by the geodesic equation, 



r = i (r-H-r) r-(M-H)- 1 -^ 



-^-H - lv(r • H • r) 



where ' = 4z and j- = ^ 



(2) 

r • V. A full solution of this 
equation (ray-tracing) is unnecessary since we may not 
only linearize the equation in H, but also linearize the so- 
lution in H by evaluating the gravitational accelerations 



on an unperturbed trajectory, i.e. the Born approxima- 
tion. A "Born trajectory" arriving at the origin at r\ = r/Q 
coming from direction n is 

r(rj) = - ^Hq • (r/ - rf) - J drf x 



H • n + i(?/ - rj) ((n • H • n) n - V(n H n 



(3) 

(r;',x') 



where x' = (770 — r/')h. The initial conditions are chosen 
such that the trajectory is light-like to the angle between 
two trajectories n and n' intersecting at the origin is 
cos -1 n • n' (both to 1st order in H). 

One may decompose the perturbation of the trajec- 
tory into a displacement along the line-of-sight (LOS), 
i.e. a time delay; and a displacement perpendicular to 
the LOS, Sr± = r (n • r) n, which describes all lensing 
effects. It is convenient to convert from displacement in 
coordinate distance to displacement in angle on the sky 
by defining A = 5 ±r / (r/o — r/) , which can be thought of as 
a 2-d vector in the tangent space of the direction sphere 
parameterized by n. 

Following [11] the displacement pattern on the direc- 
tion sphere described by A can alternately be described 
in terms of the convergence, k — — iA a :a , which is the co- 
variant divergence of A; and the rotation to — ^(A a e ab ) : b, 
which is the is the covariant divergence of A rotated by 
90° . The G and C modes are given by k and u), respec- 
tively. In constrast to the G modes, where the gravity 
wave contribution will be dwarfed by the density pertur- 
bation, the C modes will have no contribution from den- 
sity perturbations on large scales, so we are interested in 
cj and not n. Rexpressing the 2-d relation lo = ^(A a e ab )-b 
in terms of the 3-d trajectory we find 4 

w(n,rj) = -\— ^— n- (V„ x r(n,ry)) 
1 f no 

= - ^ 4/[n.(V X H).n] WiM (4) 

using the Born solution. Henceforth one can consider rj 
as being a measure of the distance to the background 
galaxies whose shear we measure. 

The contribution of the gravity waves (tensor modes) 
to H is transverse and traceless, i.e. TrH = 0, V-H = 0, 
and one can Fourier decompose the tensor constribution 
as 

^(M) = ttV / A^EHlJtM) (5) 

(2tt)2 J l= i' 



where 



Hg } (k, V ) = T (T) (k, rj) H {i) (k)e (i) (k) , (6) 



Linearized tensor (gravity wave) perturbations always lead to 
synchronous coordinates, but other coordinate choices (gauges) 
are available for vector (vorticity) and scalar (density) perturba- 
tions. 



4 One can see from this equation that a scalar (density) perturba- 
tion does not contribute to u) and hence to C modes. 
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k = |k|, wc define the transfer function such that 
T( T )(fc, 0) = 1, and i sums over the two (linear) polariza- 
tion states defined by the polarization tensors, , which 
obey Tre w = Ime (i) = 0, e (i) = e^, k ■ e (i) (k) = 0, 

Tr(e(j)(k) -e(j)(k)) = 25ij. Assuming isotropy and no 
preferred handedness we may define a power spectrum: 



HwMHfoM) = (2^) 3 P (T) (fc)^ ( 5( 3 )(k-k') . (7) 



In inflationary models with Hubble parameter H\ during 
inflation 



P (T )(fc) 
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which gives equal metric perturbation on all scales. 
Presently CMB observations limit [14] Hi < 2 x 10 14 GeV. 
An integral required to compute to is 

fi(i)(k, n, r]) = — ^ k sin 2 6 sin 2</> 
rvo 

x / dT]'e ikcos9 ^°- r i^T(k,r]) (9) 

Jq 

where the two angles, 9 and <p, are defined by n • 
k = k cos 9 and n • (k x e^^k)) • n = k sin 2 9 sin 2(f). 
Eq.s (4,5,6) becomes 

1 r 2 

w(n, r?) - —3- / d 3 k ^ % (k) (k, n, r?) (10) 
(2tt)2 7 . . 



If wc define 
then 



{r,) = J d 2 nY^ m) (h)co(h,r,) (11) 



cfin) ~E (H, m) {v)\ 2 ) 

m— — l 

= 2 J d 3 kP {T) (k) \Tf{k,ri)\ 
where, w = k(r/ — -q) and w — krj , then 



(12) 



\ 



'tt (Z + 2)! 
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^T (fc , Wo 7 w/ )M^) , (13 ) 



Like C®, but unlike £>(; iTO ), T® is invariant to a rotation 
of coordinates, and is computed most easily when k is in 
the direction of the coordinate "North Pole" . 

The quantity to gives the rotation of the apparent posi- 
tion angle (PA) wrt the coordinate grid due to the bend- 
ing of light along the LOS. This is related to the apparent 
shear of the coordinate grid, -f a b by the relation on the 



sphere (V 2 + 2) lo = — (7 ab efc c ) :ac - For this 7 a h to be in- 
dicative of the shear measured by looking at galaxy PA's 
one requires these PA's not to be preferentially aligned 
wrt to the coordinate grid. However since H at the source 
(the galaxy) is non-zero and non-isotropic, this would 
not be true if galaxy orientation were isotropically dis- 
tributed in physical space. 

Assuming physical isotropy, we must add a "metric 
shear" caused by the shearing of the coordinates wrt 
physical space, i.e. ^ a b, which is just the traceless tran- 
verse projection of — ^H. Metric shear does not cause ro- 
tation of the images, but to compare to the above we com- 
pute Alj, which is the rotation which would be caused 
by the bending of light required to produce A-f a i,: 



(V 2 + 2) Aoj 
(Vo-v) 2 ~ 



(A 7 a V) :( 
4n \ 



Vo-V 



•(VxH) T 



(14) 
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where Vj_ = V — nii-V. Correcting for the metric shear 
leads us to correct the transfer function by adding 



'I - 1 



'tt (Z + 2)! 



(Z + 2)(Z-l)y 2 (Z-2)! 
ji(w) - jt^w)) T(k,^j^) (15) 



to eq (13). Note that none of these metric corrections 
are small especially at large redshift and for the longest 
wavelengths the two terms nearly cancel. 

Physical isotropy need not be the correct approxima- 
tion as galaxies will oscillate in phase with the gravity 
waves just as a Weber bar does, and also because the 
initial galaxy shapes may have residual correlation with 
the gravity wave due to small shear they exert on the 
galaxy when it is formed. Just how the galaxy shapes 
react to these forces depends on the details of galaxy 
dynamics and formation. Here we simply ignore such ef- 
fects. Both the metric and induced shear take us away 
from the "bending of light" in flat space picture for the 
origin of shear, and while they are negligible in the case 
of scalar perturbations, they are not for the very small 
amplitude of tensor shear. 

In the observationally indicated flat FRW cosmology 
the tensor metric, H, obeys the wave equation H— V 2 H+ 
2|H = 0. With zero cosmological constant T^(k,rj) — 
3ji(krj)/(krj), and this is also very for accurate for the 
observationally indicated A, given by Q\ « 0.7. 

Figure 1 shows the spectrum of the shear produced 
by inflation. This depends on the redshift of the back- 
ground galaxies; here we have chosen a very optimistic 
source redshift of z — 3 for all background galaxies. The 
gravity wave shear 1(1 + l)Cf decreases with increasing I, 
in contrast to density perturbations which increases with 
I. The dashed curve in Figure 1 shows the results using 
only Eq. (13), i.e. ignoring the metric shear. We see that 
the metric shear correction partially cancels the signal on 



4 




FIG. 1: Expected signal in an all-sky survey in the curl mode 
in a model with Q m = 0.3, = 0.7 with (solid) and without 
(dashed) the metric shear term of Eq. (15). The Hubble rate 
during inflation which determines the amplitude of the gravity 
wave power spectrum has been set to 2 x 10 14 GeV. The 
noise estimate here assumes an all-sky survey with 1.5 x 10 10 
galaxies and the variance of the intrinsic shear equal to 0.1. 
Background galaxies are all assumed to be at fixed redshift 3. 



large scales. So a very small, virtually undetectable sig- 
nal becomes completely undetectable due to the metric 
shear. Also shown in Figure 1 is the shape noise: 

where (j 2 ) is the intrinsic rms ellipticity of the galaxies 
(iVgai in all) and / s k y is the fraction of sky covered by the 
survey. Figure 1 shows the noise for an all-sky survey 
with 100 galaxies per square arcminute, or 1.5 x 10 10 
galaxies in total, roughly the density anticipated by the 
proposed SNAP and LSST missions, though they are not 
all-sky. 

It would be wonderful if inflation-produced gravity 
waves produced a C mode of cosmic shear that could 
be detected by observing cllipticities of distant galaxies. 
Alas, the signal is too small to be detected, and the best 
hope of observing C modes remains in the polarization 
of the cosmic microwave background. 
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